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Mean-field Study of Charge, Spin, and Orbital Orderings 
in Triangular-lattice Compounds 74Ni02 (A=Na, Li, Ag) 
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Department of Applied Physics, University of Tokyo, Hongo, Bunkyo-ku, Tokyo 113-8656, Japan 

We present our theoretical results on the ground states in layered triangular-lattice com- 
pounds ylNi02 (A=Na, Li, Ag). To describe the interplay between charge, spin, orbital, and 
lattice degrees of freedom in these materials, we study a doubly-degenerate Hubbard model 
with electron-phonon couplings by the Hartree-Fock approximation combined with the adia- 
batic approximation. In a weakly-correlated region, we find a metallic state accompanied by 
\/3 X charge ordering. On the other hand, we obtain an insulating phase with spin-ferro and 
orbital- ferro ordering in a wide range from intermediate to strong correlation. These phases 
share many characteristics with the low-temperature states of AgNi02 and NaNi02, respec- 
tively. The charge-ordered metallic phase is stabilized by a compromise between Coulomb 
repulsions and effective attractive interactions originating from the breathing-type electron- 
phonon coupling as well as the Hund's-rule coupling. The spin-orbital-ordered insulating phase 
is stabilized by the cooperative effect of electron correlations and the Jahn- Teller coupling, 
while the Hund'-rule coupling also plays a role in the competition with other orbital-ordered 
phases. The results suggest a unified way of understanding a variety of low-temperature phases 
in ANi02. We also discuss a keen competition among different spin-orbital-ordered phases in 
relation to a puzzling behavior observed in LiNi02. 

KEYWORDS: multi-orbital Hubbard model, electron-phonon coupling, charge order, orbital order, metal- 
insulator transition, triangular lattice, NaNi02, LiNi02, AgNi02 



1. Introduction 

One of the most distinctive features of strongly- 
correlated electron systems is diverse cooperative phe- 
nomena. ^'^^ Correlated electron systems generally show 
intricate phase diagrams full of competing or coexisting 
states, and the phase competition often leads to exotic 
many-body phenomena. Among many aspects, there are 
two factors which promote such complexity, i.e., multi- 
ple degrees of freedom ^"^^ and geometrical frustration.^^ 
The multiple degrees of freedom are composed of charge, 
spin, and orbital of electrons. Strong electron correla- 
tions induce interplay among them, resulting in a va- 
riety of phases. In addition, the interplay often causes 
exotic response to external perturbations, such as the 
colossal magneto-resistance in perovskite manganites.^-* 
On the other hand, the geometrical frustration promotes 
the phase competition. In general, the geometrical frus- 
tration results in a huge number of low-energy degen- 
erate states, by suppressing conventional long-range or- 
ders. The degeneracy yields nontrivial phenomena such 
as complicated ordering, glassy behavior, and spin-liquid 
states. It is a long-standing problem in condensed-matter 
physics to understand a variety of phenomena which 
emerge from synergetic effects between the multiple de- 
grees of freedom and the geometrical frustration. 

The family of compounds ANi02 (A=Li, Na, Ag) is 
a typical example of such geometrically-frustrated sys- 
tems with multiple degrees of freedom. ANi02 takes two 
different lattice structures depending on the cation A, 
that is, the delafossite structure"^-* (AgNi02) and the or- 
dered rock-salt structure®) (NaNi02 and LiNi02). Both 




Fig. 1. (Color online). Schematic picture of NiOe layer in ANi02. 
NiOs octahedra share their edges with neighbors to form the tri- 
angular lattice of Ni cations. 7 denotes the Ni-Ni bond directions. 
Two eg orbitals, — and — j/^, are shown. 



structures are quasi-two-dimensional, composed of stack- 
ing of Ni, O, and A layers. The magnetic and transport 
properties are dominated by Ni cations, which are sur- 
rounded by the octahedron of oxygens. The NiOg octa- 
hedra share their edges so that the Ni sites constitute 
the frustrated triangular layers as shown in Fig. 1. Ni^+ 
cation has seven Sd electrons in the low-spin configura- 
tion: Six out of seven fully occupy the lower t2g levels 
and the remaining one electron enters in the higher Cg 
levels. Hence the doubly-degenerate Cg orbital degree of 
freedom is active in these systems. 

ANi02 shows a variety of behaviors depending on the 
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cation A in spite of the similarity in lattice and electronic 
structure. NaNi02 is a Mott insulator with a gap of 
0.24cV.^^ This compound shows a first-order structural 
transition accompanied by cooperative Jahn- Teller dis- 
tortions at and a second-order magnetic tran- 
sition at 20K.i°"i2) At the lowest temperature (T), the 
system exhibits orbital-ferro and A-type antiferro-spin 
order (antiferromagnctic stacking of spin-fcrro ordered 
layers). ^"^^ LiNi02 is also a Mott insulator with a gap 
of 0.2eV,^^' however, it shows no clear phase transition 
down to 1.4K in contrast to NaNi02.^^-' There is strong 
sample dependence in T dependence of the magnetic sus- 
ceptibility; the field-cool and zero-field-cool bifurcation 
appears in different ways depending on samples. This 
sample dependence strongly suggests the relevance of 
extrinsic disorder. The ground state as well as the 
finite-T properties remain controversial. ^^"^^^ In addi- 
tion to these insulating materials, recently, a new metal- 
lic compormd AgNi02 was synthesized.^' It shows a 
structural transition associated with a -s/S x \/3 charge 
ordering at 365K and antiferromagnctic transition at 
20K.^^'^^) The system remains metallic down to the low- 
est T. It was claimed that in the low-T phase the system 
separates into rather localized spins at 1/3 Ni sites and 
itinerant electrons at the remaining sites. Magnetic 
properties at low T were analyzed by considering the 
competing nearest- and second-neighbor exchange cou- 
plings between localized spins. ^^'^^^ 

So far, ANi02 has been theoretically studied by the 
first-principle calculations and the strong-coupling anal- 
yses. The orbital and spin ordering observed in NaNi02 
was reproduced by the LSDA-I-J7 first-principle calcula- 
tions. '^°''^^'' The metallicity as well as the charge ordering 
in AgNi02 was also reproduced by the first-principle cal- 
culations.^^' On the other hand, effective models in the 
limit of strong electron correlation, the so-called Kugel- 
Khomskii models, were studied to understand the orbital 
and spin ordering in NaNi02 and the peculiar disordered 
state in LiNi02.^°'^^'^^^ For AgNi02, recently, the mag- 
netic phase diagram was investigated by the classical spin 
model with ignoring the itinerant electrons.'^*-' 

Despite the extensive studies so far, comprehensive 
understanding of the ground states of ANi02 has not 
been reached yet. Although the low-T states of NaNi02 
and AgNi02 are reproduced by the first-principle calcu- 
lations, the mechanism of stabilizing these states is not 
fully clarified. In addition, the effective-model approach 
does not fully succeed in reproducing the ground states 
of NaNi02 and LiNi02 in an unified way. A possible way 
to explore the comprehensive understanding is to inves- 
tigate a model in a wide region of interaction parameters 
systematically beyond the strong-coupling approach. In 
fact, in both NaNi02 and LiNi02, the Mott gap is not 
large; the gap is comparable to the transfer integrals. Fur- 
thermore, the newly-synthesized AgNi02 shows metallic 
behavior. These facts suggest an importance of charge 
fluctuations in weakly- or intermediately-correlated re- 
gions. Electron-phonon couplings may be another essen- 
tial factors, which have not been considered seriously in 
spite of the experimental facts that the structural tran- 
sitions are observed in NaNi02 and AgNi02. 



In this study, aiming at a unified picture of these 
compounds ANi02, we investigate a multi-orbital Hub- 
bard model with electron-phonon couplings on a two- 
dimensional triangular lattice. Our purpose is to eluci- 
date the microscopic mechanism for the variety of phases 
in these compounds ANi02. In particular, we focus on 
electron-phonon couplings and charge degrees of free- 
dom, both of which have not been carefully examined 
in the previous studies. We clarify that these two ele- 
ments play an important role in the phase competition 
in these complicated systems. We obtain a y/3 x -y/S 
charge-ordered metallic (COM) phase in the weakly- 
correlated region and an insulating phase with ferro- 
type spin and orbital ordering in the intermediatcly- 
to strongly-correlated region. These two phases repro- 
duce many aspects of the low-T states in AgNi02 and 
NaNi02, respectively. We discuss the peculiar disordered 
state in LiNi02 in relation with a keen phase competition 
in the obtained phase diagram. 

The organization of this paper is as follows. In Sec. 2, 
we describe our model and method. We introduce the 
Hamiltonian term by term, and present the approxima- 
tions adopted in the calculations. We present our results 
in Sec. 3. We discuss the parameter region and the mech- 
anism to stabilize the a/3 x V3 COM phase and the 
spin-orbital-ordered insulating phase. Finally, Sec. 4 is 
devoted to summary. 

2. Model and Method 

2.1 Model 

In the present study, to elucidate the phase competi- 
tion in ANi02, we investigate the ground state of the 
multi-orbital Hubbard model with electron-phonon cou- 
plings. Among the five id orbitals, we consider only the 
twofold degenerate eg orbitals, by taking account of the 
low-spin state of Ni^"*" cations. Our Hamiltonian is writ- 
ten as 

"H = Ti-kin + ^int + ^ol-ph + ^ph, (1) 

where "Hkin, "Hint, ^ci-ph, and Hph represent the ki- 
netic term of electrons, the electron-electron interactions, 
the electron-phonon couplings, and the elastic term of 
phonons, respectively. We describe the detailed forms of 
each term in the following. 

2.1.1 Kinetic term 

Due to the spatial anisotropy of the e^-orbital wave 
functions, transfer integrals between Ni sites depend on 
the bond direction as well as orbital types (see Fig. 1). 
The kinetic term in eq. (1) is written as 

^kin = - E E E {4c..C,p. + H.C.) . (2) 
(tj) a 

Here, i and j denote the site indices, a and /3 represent 
the orbital indices with a = a(6) corresponding to the 
3^2 _ ^2 ^^2 _ y2-^ orbital, a is the spin, and 7,,, denotes 
the direction of bond between the site i and j, as shown in 
Fig. 1. The sum over (ij) is taken for the nearest-neighbor 
sites on the triangular lattice. The transfer integrals are 
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given by the following matrices; 



t 
t' 



t2 t3 



t2 -ts 



(3) 

for the two bases of iz^—r^ and x^—y^ orbitals. From the 
symmetry of orbitals, we obtain the following relations: 
t2 = t/A + 3t'/A, ts = V^{t-t')/A, and = 3i/4 + t'/4, 
with two independent parameters, t and t' . We set t — 1 
as an energy scale. The value of t' depends on both d- 
d direct transfer integrals and d-p-d indirect ones in a 
complicated manner. '^^''^^^ In the following, we show the 
results for t' = —1 hy noting that the orbital overlaps 
between atomic orbitals at the neighboring sites lead to 
t' /t — 1 when one consider both contributions. An ex- 
tended study in wider range of f /t for an effective model 
without phonon is found in Ref. 33. The choice of t and 
t' gives the non-interacting bandwidth 8i. 

2.1.2 Electron- electron interactions 

Next we introduce the electron-electron interaction 
term "Hint in cq. (1). We consider only the on-site 
Coulomb interactions. For the doubly-degenerate eg or- 
bital system, "Hint is written as 



Hint ^Hu+ Hu' + Hj^ + Hj' , 



where 



i a 



'T-iacr'T-ibo.' ) 



(4) 

(5) 
(6) 

(7) 

(8) 



Here Uiaa = c[aa^iaa, U and U' denote the intra- and 
inter-orbital Coulomb repulsions, and Jh and denote 
the exchange interaction and the pair hopping, respec- 
tively. "Hjjj and Hj^ are called the Hund's-rule cou- 
plings. Hereafter we assume the relations U' = U — 2 Jh 
and Jh = Jh to retain the rotational symmetry of the 
Coulomb interaction. 

2.1.S Electron-phonon couplings 

As to the electron-phonon couplings, we consider two 
relevant distortions of NiOg octahcdra in ^Ni02 , namely, 
the Aig breathing mode and the Eg Jahn- Teller modes. 
Hoi-ph in eq. (1) is given by the sum of these two contri- 
butions as 



i-cl-ph ^ "-cl-ph- (9) 

The Aig mode corresponds to the isotropic expan- 
sion (contraction) of NiOe octahedron [Fig. 2(a)], which 
lowers (raises) two Cg energy levels without lifting their 
degeneracy. Namely, the Aig mode couples to the local 
charge on each Ni site, written in the form 



K 



el-ph 



7br E ^br,i (n-ia + ni& - 1) , (10) 



(a) 



(b) 



I 








Fig. 2. Schematic pictures of the displacement of oxygens (open 
circles) around a Ni cation at the origin in (a) Aig mode, (b) 
Eg^u mode, and (c) -Eg,ii mode. 



where Uia = X^ct^mo-j ^br,i is the amplitude of the Aig 
lattice distortion, and 7br > is the corresponding cou- 
pling constant. A positive (negative) Xbr.i corresponds to 
an expansion (contraction). 

The Eg mode has two components, Eg^u and Eg,^, as 
shown in Figs. 2(b) and 2(c), respectively. The Eg^u mode 
corresponds to the z-axis elongation of NiOe octahedron, 
which splits the energy levels of x'^ — and — r^ 
orbitals, while the Eg^y mode causes a mixing of the two 
orbitals: The coupling to the Eg modes is written as 



n. 



JT 
cl-ph 



-7JT E{^JT,i [riiaa - nibcr) 



+ a;jT,j(cL^Q6<, -I- 4^C,a<T)}, (11) 

where xjT.i and xjT,i are the amplitudes of Eg_u and 
Eg^^ modes, respectively, and 7jt represents the common 
coupling constant. 

2.1.4 Phonon term 

The phonon term Hph in eq. (1) consists of the on-site 
term {W^^^''') and the inter-site term (Hp^"'') as. 



nj -T/clastic 

n-ph — 'T-ph 



^coop 

^ph • 



(12) 



Each term is given by the sum of contributions from the 
Aig and Eg mode phonons. The first term T^pj^*^*"^ is the 
elastic energy of lattice distortions, which is given by the 
sum of the following two terms; 



-T/clastic,br 1 X ^ 2 



clastic, JT 

ph 



-E 
2 ^ 



'JT,j: 



'-JT 



(13) 



(14) 



The elastic constants of Aig mode and Eg modes are 
taken as unity without losing generality, by normaliz- 
ing the amplitudes of lattice distortions (xbr,i, a^jT.i, and 
xjT.i) and the coupling constants (7br and 7jt)- 

The second term V.'^^ describes the cooperative cou- 
plings of lattice distortions. The precise form of the cou- 
plings is, in principle, determined by the phonon disper- 
sion in the materials, but here, for simplicity, we con- 
sider the nearest-neighbor couplings only. Then 'H'^^ is 
defined by the sum of two terms; 



n 



coop.br 



ph 



= Abr E ^^'^ 



(15) 
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coop. JT 

ph 



\JT 7 ,a;jT,ia;jTj, 



(16) 



with the couphng constants Abr, Ajt, and Ajt- Although 
the values of Ajt and Ajt are generally different, we take 
Ajt = Ajt for simplicity. It is reasonable to assume Abr 
to be positive, i.e., the "antiferro" -type coupling, since 
the ^ig-type expansion (contraction) of NiOg octahe- 
dron tends to shrink (expand) the neighboring octahedra 
due to the edge-sharing network of octahedra. On the 
other hand, the tendency is opposite for the Eg modes; 
the Jahn- Teller distortion of an octahedra induces the 
same distortion in the neighboring octahedra. Hence we 
consider the "ferro"-type coupling Ajt > in the follow- 
ing study. In the following calculations, the "antiferro" - 
type Aig coupling tends to stabilize a charge ordering by 
differentiating charge density between the neighboring 
sites, while the "ferro"-type Jahn- Teller coupling favors 
"ferro"-type orbital ordering. 

2.2 Method 

In order to study the ground state of the model (1) 
in a wide range of parameters, we adopt the Hartree- 
Fock approximation to decouple the electron-electron in- 
teractions, and the adiabatic approximation to treat the 
electron-phonon couplings. Within the Hartree-Fock ap- 
proximation, the two-body interaction terms in Hint are 
decoupled by introducing mean fields, {c\aa^ia' a') ■ The 
amplitudes of lattice distortions are determined by the 
adiabatic approximation. Within this approximation, the 
equilibrium values of Xbr.i, a;jT,i, and a;jT,i are deter- 
mined by using the Hellmann-Feynman theorem as 

dH \ / dU \ / dU \ „ ^^^^ 



dXhr. 



dx 



31A 



dx 



JT,i 



0. 



These relations lead to the set of equations in the form 

Xhr.i + Abr^ Xhr,j =7br ((?^^a) + {nib) - 1) , (18) 
j 

XjT.i - AjT^ a;jTj =7.IT ^ {{Uiaa) - {niba)) , (19) 

SjT.i - -^JT^ XjT^j ^liT^{{c\acrCiba) + (cLc^aa)), 
3 

(20) 

where the sum is taken over the nearest neighbors of 
the site i. 

We determine the mean fields and the lattice dis- 
tortions in a self-consistent way. For a given set 
of {(cL<TCjaV'):2;br,j,a;jT,i,:^jT,J, we diagonalize the 
Hamiltonian under the Hartree-Fock approximation and 
obtain one-particle eigenenergies and eigenstates, which 
are used to calculate a new set of {0iaa'^ia' a')} ■ These 
new mean fields are substituted in Eqs. (18)-(20) to de- 
termine the new set of {xbr.i, a;jT,i, a;jT,i}- These pro- 
cedures are repeated until the convergence is reached 
within the precision less than 10~* for all the variables. 

In the following calculations, we take the unit cell 
which includes six Ni^+ sites in the triangular lattice 



as shown in Fig. 3. To incorporate different orbital or- 
derings depending on the way of taking the 7 direction 
in Fig. 1, we consider two different ways of embedding 
the unit cell as shown in Fig. 3. These unit cells accom- 
modate a a/S X \/3 charge ordering and a two-sublattice 
ordering such as a stripe-type antiferromagnetic state. 
Note that the ordering patterns observed in NaNi02 and 
AgNi02 are both included by taking the unit cells. For 
the initial state in the iteration, we consider more than 30 
different states with different symmetry in spin, orbital, 
and charge sectors, which are relevant in the parameter 
space we study. The different initial configurations are 
adopted for each type of the unit cells in Fig. 3. 

The integration over the wave number in the calcula- 
tions of the mean fields is replaced by the sum over 24x 24 
grids in the Brillouin zone for the supercell. Hereafter, 
we focus on the quarter-filling case, i.e., one electron per 
site on average, corresponding to one Cg electron in the 
low-spin state of Ni3+. 




'y = 3 



unit cell 



Fig. 3. Two different ways of taking the unit cell with six Ni sites 
on the triangular lattice used in the calculations. The axis 7 for 
the orbitals is shown (see also Fig. 1). 



3. Results and Discussion 

In this section, we show the results obtained for the 
ground state of the model (1). In particular, wc focus on 
the roles of U . J-a/U, and 7jt, and discuss how these 
parameters affect the ground state. We take 7br = 1.6 
and Abr = Ajt = 0.05, which result in reasonable energy 
gain in forming CO, as we will see below. 

As a result, we find that the tendency to charge or- 
dering becomes pronounced by a compromise among 
U, U' , Jh, and 7br. On-site repulsions U and U' sup- 
press charge disproportionation, while the Hund's-rule 
coupling Jh as well as the breathing-type coupling 7br 
works as an inter-orbital effective attractive interaction 
and promotes charge disproportionation. On the other 
hand, for larger U and U' , the system becomes insulat- 
ing, and the Jahn- Teller coupling 7jt becomes important 
and enhances the tendency to orbital ordering, concomi- 
tant with magnetic ordering. In order to characterize the 
orbital-ordered phases, we introduce the pseudospin op- 
erators in the orbital sector, defined as 



(21) 



where cr^^ denotes the pauli matrix. For example, Tz-OF 
means a ferro-type order of t^^ , which is the z component 
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of Tj. 

The main result is summarized as the phase diagrams 
shown in Fig. 4. In the following, we will focus on the 
\/3 X "v/S-type charger-ordered metallic (COM) phase 
found in the weak-coupling region [Fig. 4(a)], and the 
orbital-ferro spin-ferro insulating (r^-OF SF I) phase ob- 
tained in the intermediately- to strongly-correlated re- 
gion [Figs. 4(b) and 4(c)]. These two phases are candi- 
dates for the low-T states of AgNi02 and NaNi02; re- 
spectively. We will identify the parameter range for these 
phases, and discuss the origin and the nature of them in 
the following. 

3.1 Weakly- correlated region 
3.1.1 Phase diagram 

A tendency toward charge ordering is widely observed 
in a weakly-correlated region. Figure 4(a) shows the 
phase diagram at U = 4. Among the competing phases, 
a COM phase with \/3 x VS-type charge ordering is sta- 
bilized in a wide range of 7jt at Jn/U ~ 0.25. The 
X >/3 charge ordering is a three-sublattice order, in 
which one is charge rich (the density is almost two) and 
the other two charge poor (the density is almost 0.5 per 
site). This \/3 x -y/S COM state, shown in Fig. 5(a), is 
remarkable, since it has the same charge ordering pat- 
tern as the low-T state of AgNi02.^'''^'^^ The spin state 
is also interesting; large moments appear at charge-rich 
sites (5 ~ 0.6), while the moments are suppressed at 
charge-poor sites {S ~ 0.05). A similar differentiation 
was proposed for the low-T state of AgNi02,^^' al- 
though the calculated spin pattern does not fully agree 
with the experimental result. We note that importance 
of interlayer coupling is experimentally suggested for the 
magnetic ordering, ^^"^^^ which is not taken into account 
in our model. We will discuss the nature of this COM 
phase in Sec. 3.1.3 in detail. 

Around Jb_/U = 0.24 and 7jt ~ 1-1, we find another 
COM phase, i.e., the sixfold COM phase. Although this 
phase is stabilized in a narrow region in the phase di- 
agram, it remains to be metastable in a wide parame- 
ter range, with a slightly higher energy than the ground 
state, as we will discuss in the next section 3.1.2. This 
phase has the consistent ordering structure with the low- 
T state of AgNi02 in terms of both charge and spin, as 
shown in Fig. 5(c). Strictly speaking, the charge pat- 
tern of this phase has lower symmetry compared with 
AgNi02, due to the superposition of stripe- type charge 
modulation onto the -s/S x a/3 charge ordering. How- 
ever, the magnitude of this modulation is very small: 
For charge rich sites (sites 1 and 4 in Fig. 3), we have 
(ni) = 1.5040 and (77,4) = 1.5027, while for charge 
poor sites (sites 2, 3, 5, and 6 in Fig. 3), we have 
(na) = (ns) = 0.7479 and (na) = (jie) = 0.7487, at 
Jn/U = 0.24 and 7jt = 1.1. The modulation gives very 
small charge disproportionation of the order of ^ 0.001 
within charge rich and poor sites. [The small dispropor- 
tionations are exaggerated in the schematic picture in 
Fig. 5(c).] 

In addition to the two COM states, we obtain a va- 
riety of ordered phases in the weakly-correlated region. 
Among them, we focus on two phases which compete 
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Fig. 4. (Color online). Ground-state phase diagrams for the 
model (1) in the plane of J^/U and 7jt at (a) (7 = 4, (b) 
{/ = 8, and (c) U = 20. We take t' = -1, 7br = 1.6, and 
-^br = -^.TT = 0.05. The ordering patterns of each phase arc 
shown in Fig. 5. 



with the COM states; the sixfold charge-ordered insulat- 
ing (sixfold COI) phase in the large Jn/t/ region, and 
the spin-ferro metallic phase with a weak charge order- 
ing (SF-COM) stabiHzed for smaller Ju/U [Fig. 4(a)]. 
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(a) 73 X 73 COM (b) OF SF I 




(e) OFSAFM(I) (f) sixfold COI 




(g) OAF-COI 1 (h) OAF-COI 2 




Fig. 5. Schematic pictures of charge, spin, and orbital ordering 
patterns in the six-site unit cell for the phases in Fig. 4. The 
size of the circles schematically indicates the charge density at 
each Ni site. The length of the arrows denotes the magnitude 
of the spin moment at each site. The open circles denote the 
orbital para state, while the filled or shaded circles show orbitally 
polarized states. The different patterns denote different orbitally 
polarized states. See the text for details of (a)-(d). Both tz-OF 
SF I and r^-OF SF are represented by Fig. 5(b) because these 
phases have the same symmetry (OF SF). Both OF SAF I and 
Tx-OF SAF M (1) are represented by Fig. 5(e). Typical values 
of of the charge disproportionation and orbital polarization in 
(e)-(i) are as follows, (e) OF SAF 1: (t^) ~ (t^) ~ 0.6 at all 
sites. r^-OF SAF M: {t^) ~ 0.2 and (r^) ~ -0.06 at all sites, (f): 
(n) ~ 3.0 and (n) ~ 1.7 at the two charge-rich sites and (n) ~ 0.3 
at the other charge-poor sites, (g): (tx) ~ 0.2, (tz) ~ 0.7, and 
(n) ~ 1.2 at the charge rich sites; (tx) ~ 0.1, (r^) ~ —0.2, and 
(n) ~ 0.8 at the charge poor sites, (h): (tx) ~ 1.4, (t^) ~ 0.6, 
and (n) ~ 1.7 at the charge rich sites and (tx) ~ 0.2, (r^) ~ 0.0, 
and (n) ~ 0.3 at the charge poor sites, (i): (n) ~ 1.02 at the 
charge rich sites (1 and 4 in Fig. 3), (n) ~ 0.98 at the charge 
poor sites (3 and 6), and (n) ~ 1.00 at the other charge-neutral 
sites (2 and 5). 

We argue the stability of the -\/3 x \/3 COM as weU as 
the sixfold COM in comparison with the two competing 
phases in the next section. 

3.1.2 Stability of the x VS COM phase 

In order to clarify the competition among the ^/3x ^3 
COM, sixfold COM, sixfold COI, and SF-COM phases, 
we investigate the internal energy in detail by comparing 
the contributions from different terms in the Hamilto- 
nian; E]^in = (Hkin), Eu = CHu), Eu' = CHu'), Ej^ = 
(H,7h +nj^), iSbr ^ and 
the total energy Etot = (H) . We show the comparison as 



a function of U at Jn/U = 0.27 and 7jt = 0.5 in Fig. 6. 

For U < 3.9, the sixfold COI state has the lowest en- 
ergy. As shown in Fig. 5(f), this phase has a polaronic 
nature, namely, one site is almost fully occupied (the 
local density is almost 4), and another one site accom- 
modates almost two electrons. At the latter site, spins 
of two electrons are aligned parallel by the Hund's-rule 
coupling. This phase is stabilized in a region where the 
repulsive Coulomb interactions are compensated by effec- 
tive attractive interactions originating in the breathing- 
type electron-phonon coupling as well as the inter-orbital 
Hund's-rule coupling. In fact, it is clearly observed in 
Figs. 6(d) and 6(e) that the energy gain in Ej-^ and £'br 
contributes to the stabilization of the sixfold COI phase. 

On the other hand, for U > 4.9, the SF-COM state 
is most stabilized. As schematically shown in Fig. 5(d), 
the charge ordering in this phase is a stripe type, but the 
charge disproportionation is very small ((n) ^ 1.03—1.09 
at charge rich sites, while (n) ~ 0.91 — 0.97 at charge poor 
sites): the main feature is the spin ferromagnetic order- 
ing. The origin of this phase can be attributed to the 
Stoner mechanism. As shown in the inset of Fig. 7, 
the non-interacting Fermi level is located in the vicinity 
of the steep peak of the density of states (DOS). Conse- 
quently, a ferromagnetic instability is caused at a rela- 
tively small U ~ 3.4, according to the Stoner's criterion. 
The characteristics of Stoner ferromagnet are observed in 
the energy comparison in Fig. 6(c); Eu becomes smallest 
among the competing phases. 

The VSx -\/3 COM state intervenes these two, and has 
the lowest energy for 4.0 < 4.8. In the same parame- 
ter range, the sixfold COM state appears as a metastable 
state and stays very close to the ground state, as shown 
in the inset of Fig. 6(a) (the energy difference is less than 
0.02). These two COM states are stabilized by a compro- 
mise between the different stabilization mechanisms for 
the sixfold COI and the SF-COM phases. According to 
Fig. 6, the COM phases have higher (lower) Ej-^ and 
£^br, while they have lower (higher) Ejj and Ei/i, com- 
pared with the sixfold COI phase (the SF-COM phase). 
Namely, the two COM phases are stabilized in a subtle 
balance between repulsive Coulomb interactions and ef- 
fective attractive interactions due to the Hund's-rule cou- 
pling and the breathing-type electron-phonon coupling. 

Since the COM phases are stabilized in a delicate com- 
promise, it is important to consider their stability against 
the elements which are ignored in our current analysis, 
such as fluctuations beyond the mean-field level and the 
long-range part of electron interactions. First, we note 
that the Stoner ferromagnetism is fragile when consider- 
ing the electron correlation effect beyond the mean-field 
approximation. •^^^ Therefore, the COM phases are ex- 
pected to extend to larger U or smaller J-a/U . Second, 
the amplitudes of breathing-type distortions are fairly 
large in the sixfold COI phase compared to those in the 
other phases. Hence this phase will be suppressed by con- 
sidering more realistic contributions from phonons, e.g., 
anharmonic terms of phonons. This may give a chance 
for the COM phases to become wider also in smaller 
U or larger Jti/U regions. Finally, the long-range part 
of electron interactions generally works in favor of the 



J. Phys. Soc. Jpn. 



Full Paper 



Author Name 7 



(a) sixfold COI -^73 X 73 COM-»-^ SF-COM 




3 4 5 6 

U 



Fig. 6. (Color online). Energy comparisons among the X 
COM (circle), sixfold COM (triangle), sixfold COI (square), and 
SF-COM states (diamond): U dependences of (a) the total en- 
ergy, and the contribution from (b) kinetic term, (c) Coulomb 
repulsions [closed (open) symbols denote Eu (Eur)], (d) Hund'- 
rule coupling, and (e) breathing-type electron-phonon coupling. 
The inset of (a) shows the energy differences between the \/3x \/3 
COM phase and the other competing phases. The parame- 
ters are t' = —1, Jyi/U = 0.27, 7br = 1-6, 7jt = 0.5, and 
Abr = AjT = 0.05. 



charge ordering, in particular, the -y/S x type and the 
sixfold COM, as is discussed in several transition metal 
compounds and organic materials.'^''' Therefore, we ex- 
pect that the COM phases relevant to AgNi02 become 
more stable in a wider parameter range when extend- 
ing the analyses beyond the present model and method. 
Although the COM phases are robust in this parameter 
region, the energy difference between the sixfold COM 
phase and the -s/S x VS COM phase is very small, im- 
plying that the magnetic ordering pattern might be af- 
fected by small perturbations, such as inter-layer cou- 
pling. More accurate studies beyond the mean-field ap- 
proximation are necessary for fully determining the spin 
pattern of the ground state. 



3.1.3 Nature of the VS x COM phase 

Reflecting the subtle balance between the attractive 
and repulsive interactions, the -s/S x ^/S COM phase 
shows peculiar electronic properties. The density of 
states (DOS) in the y/3 x -\/3 COM phase is shown 
in Fig. 7. The site- and spin-resolved DOS indicates 
that the system exhibits a half-metallic nature: Up-spin 
electrons are localized at charge-rich sites, showing a 
gap at the Fermi level, on the other hand, down-spin 
electrons remain conductive, with a finite DOS at the 
Fermi level. Electron correlations dominantly affect up- 
spin electrons; down-spin conductive electrons preserve 
the non-interacting band structure. For comparison, we 
show DOS for the non-interacting case in the inset of 
Fig. 7, which is quite similar to DOS of down-spin con- 
ductive electrons. DOS in the VSx ^/S COM phase repro- 
duces several aspects of the result obtained by the first- 
principle band calculation:^^) The electrons at charge 
rich site tend to localize and the electronic structure at 
the other two charge poor sites resembles each other. 

This peculiar electronic state can be attributed to the 
delicate balance between the attractive and repulsive in- 
teractions. Wc plot the effective one-body potential in 
Fig. 8, which is defined as the sum of the terms in "Hint 
and Hoi-ph, which couple to the density operator rii at 
each site under the Hartree-Fock approximation. Fig- 
ure 8 shows that the charge-rich (-poor) sites bear at- 
tractive (repulsive) potentials for up-spin electrons. In 
contrast, the cancellation between the breathing-type 
electron-phonon coupling and the repulsive Coulomb in- 
teractions leads to an almost fiat potential for down-spin 
electrons. Consequently, these interactions only work as 
a shift of chemical potential, and the down-spin electrons 
retain the non-interacting band structure. 

To conclude the discussions for the weakly-correlated 
region, the -y/S x \/3 COM phase is stabilized by a com- 
promise between repulsive Coulomb interactions and at- 
tractive interactions originating from the breathing-type 
electron-phonon coupling as well as the Hund's-rulc cou- 
pling. This phase shows a half-metallic behavior with 
large magnetic moments almost localized at charge-rich 
sites and conduction electrons moving almost freely in 
the entire lattice. We note that this phase is distinguished 
from the so-called pinball liquid state, in which the elec- 
trons at the charge-rich sites exclude the conduction elec- 
trons as hard core potentials and confine them to the 
honeycomb network of charge-poor sites, as discussed in 
a spinless tight-binding model with intersite Coulomb re- 
pulsion on the triangular lattice. "^^^ 

3.2 Intermediately- to strongly- correlated region 

Next, let us consider the intermediately- to strongly- 
correlated region. Representative phase diagrams are 
shown in Figs. 4(b) and 4(c). 

3.2.1 Phase diagram 

We show the ground-state phase diagram at ?7 = 8 in 
Fig. 4(b). Note that the value of U is comparable with 
the non- interacting bandwidth., i.e., the system is in the 
intermediately-correlated region. In this region, we find 
three dominant orbital-ordered phases; Tj,-ordered spin- 
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Fig. 7. (Color online). DOS per site for the ^3 X COM phase. 
The Fermi level is set to be the origin of energy. As to the \/3x \/3 
COM phase, the up-spin (down-spin) DOS is drawn on the upper 
(lower) side. The parameters are chosen as [/ = 4, Jfi/U = 0.27, 
7br = 1-6, 7JT = 0.5, and Abr = Ajt = 0.05. The inset shows 
DOS for the non-interacting case. 



understand the origin of these phases from the strong- 
coupUng analysis, i.e., by starting from the Mott insulat- 
ing state at [/ = oo. 

In fact, the Tj,-OF SF I state was obtained for an effec- 
tive spin-orbital model in the strong-coupling limit (F3 
phase of Fig. 8 in Ref. 33). Our result is consistent with 
the previous study. Meanwhile, the competition between 
the T^-OF SAF I phase and the r^-OF SF I phase is 
obtained for the first time by explicitly taking account 
of electron-phonon couplings. In the following, we will 
consider the mechanism of stabilization of these phases 
through the detailed study of energetics. 



3.2.2 Stability of the Tz-OF SF insulating phase 

In order to understand the stability condition, it is 
instructive to rewrite H^iipi^, Hj^, and "Hj^ by using the 
pseudospin operators in eq. (21) as 



njJT 
''•el-ph 



{riia + nib)}. 



(22) 



(23) 




site 



Fig. 8. (Color online), (a) Effective one-body potentials at each 
site for up-spin and down-spin electrons (see the text for details) . 
The horizontal axis denotes the site indices in the unit cell, shown 
in Fig. 3. The charge-rich (-poor) sites are the site 1 and 4 (2, 
3, 5, and 6). The parameters are chosen as !7 = 4, J^/U = 0.27, 
7br = 1-6, 7JT = 0.5, and A^r = A^r = 0.05, consistent with 
Fig. 7. 



ferro insulator (ry-OF SF I), r^^-ordered spin-antiferro 
insulator [tx-OY SAF I), and r^-ordered spin- ferro insu- 
lator (tz-OF SF I). The ordering pattern of Ty- or Tz-OY 
SF I (r^-OF SAF I) is schematically shown in Fig. 5(b) 
[Fig. 5(e)]. For the small 7jt and Jn/U region, Ty-OF 
SF I is stabilized, while it is replaced by Tx-OF SAF I for 
larger 7jt or by a charge-ordered state for larger Jh/U. 
Meanwhile, when both 7jt and Jn/U become large, r^- 
OF SF I is stabilized. Among these phases, the r^-OF 
SF I phase deserves attention, since the spin and orbital 
pattern of this phase is consistent with the low-T phase 
of NaNiOa. 

Remarkably, the three OF phases remain stable in a 
wide range of U toward the strongly-correlated regime. 
Figure 4(c) shows an example of the phase diagram at 
large U. The result indicates that the three phases remain 
in similar parameter regions of 7jt and Jh/U compared 
to the intermediate- f7 case. This implies a possibility to 



These equations clearly show that the Jahn- Teller cou- 
pling stabilizes the Tx and Tz orbital orderings, while the 
Hund's-rule coupling destabilizes the Tx orbital ordering. 
In Fig. 9, we compare the energy contributions including 
these terms, Ejt ^ {H^^^ + H'^-^'''^'^ + and 
Ej^ , together with other relevant energy contributions, 
for the three orbital-ordered phases. 

Figure 9(b) shows that the stability of the t^-OF SF I 
phase in the small 7jt region is attributed to the energy 
gain in the kinetic energy ii^kin- This is consistent with 
the result of strong-coupling analysis, where the kinetic 
energy gain through the spin-orbital superexchange in- 
teractions is claimed to be the origin of this phase. 
It is also observed that all energy contributions in this 
phase are fairly insensitive to 7jt, as shown in Fig. 9. 
This is expected from the absence of coupling between 
Ty and Jahn- Teller distortions, as is clear from eq. (22). 

On the other hand, the r^^^-OF SAF I and r^-OF SF 
I states lower their energy through the coupling to the 
Jahn- Teller distortions [Fig. 9(e)], as expected from eq. 
(22); thus they replace the Tj,-0F SF I phase and become 
the ground state for larger 7jt, as shown in Fig. 9(a). The 
Hund's-rule coupling plays an important role in the rel- 
ative stability between the ra;-OF SAF I and the r^-OF 
SF I phases. As is evident from eq. (23), the Hund's-rule 
coupling affects only the Tx ordering, and destabilizes it 
[Fig. 9(d)]. Furthermore, the r^-OF SF I phase is stabi- 
lized by the kinetic energy gain from the interorbital hop- 
pings, compared with the Tx-OF SAF I phase as shown in 
Fig. 9(b). In fact, according to the second-order pertur- 
bation from the strong coupling limit U — >■ oo, the r^-OF 
SF I phase has lower energy than Tj,-OF SAF I phase for 
Jh/U > ^ 0.27. The phase boundary between the 

two phases in Fig. 4(c) is roughly located around this 
critical value, which indicates that the phase competi- 
tion is essentially understood from the strong coupling 
picture. 
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Fig. 9. (Color online). 7,jt dependences of (a) the total energy, 
and the eontribution from (b) kinetic term, (c) Coulomb re- 
pulsions [closed (open) symbols denote Eu [Eui)\, (d) Hund's- 
rule coupling, and (e) Jahn- Teller coupling. Comparison is made 
for the T^-OF SF I (circle), r^-OF SAF I (square), and r^-OF 
SF I states (triangle). The parameters are t' = —1, U = 20, 
Jy^/U = 0.25, 7bi. = 1.6, and Abr = Ajt = 0.05. 



3.2.3 Nature of the t^-OF SF insulating phase 

Due to the strong electron repulsion, an excitation gap 
opens at the Fermi level in the Tz-OF SF I phase. Hence, 
this state is insulating, consistent with the low-T insu- 
lating phase of NaNi02. Figure 10 shows DOS in the 
r^-OF SF I phase. DOS is composed of four sectors (two 
in each spin component), and the total weight of each 
sector is equal to one electron per site, as indicated in 
the integrated DOS in the figure. This structure can be 
qualitatively understood by considering the atomic limit 
with ignoring the hoppings, t and t' . Let us assume the 
perfect r^-OF SF order, i.e., ('T.i,3z2_^2^) = 1 and let 
other mean fields to be zero. Then the excitation ener- 
gies in the atomic limit are estimated as 

E,2_y2^^^U' - Jh~2E:^^, (24) 

i?,2_,2,^ = C/'-2i?*T, (25) 

£^3.= -r^; = C/ + 2£;j*T, (26) 

where Ea^y signifies the energy necessary to add one elec- 
tron with orbital a and spin a to the (30^ — r^, f) ground 
state. i?jrp is the energy gain from the Jahn- Teller dis- 



tortion per site, estimated as 

Substituting the parameters used in Fig. 10 (17 = 20, 
U' = 10, Jh = 5, 7jt = 1-0, and Ajt = 0.05) into these 
equations, we obtain E^2_y2^^ ~ 6.4, E^2_y2^^ ~ 11, and 
i?3j,2_^2^ ~ 19. These values well correspond to the mean 
energy of each sector of DOS measured from that for the 
lowest one in Fig. 10. 

Eq. (24) gives an estimate of the energy gap A in the 
atomic limit. This atomic value is reduced for finite t and 
t' , since each atomic level is broadened by a renormal- 
izcd bandwidth W: The estimate of energy gap is mod- 
ified as A ~ E^2_y2^^ — W . We use this simple estimate 
with replacing W by the bare bandwidth W = St for a 
comparison with the Hartree-Fock solutions. As shown 
in Fig. 11, our simple estimate from the atomic limit is 
qualitatively consistent with Hartrcc-Fock results. This 
fact also supports that the electronic spectrum of the 
r^-OF SF I phase is adiabatically continued from the 
strong-coupling limit. 

To conclude this part, three orbital-ordered insulat- 
ing phases appear dominantly in the intermediately- to 
strongly-correlated region. Among them, the T2-OF SF 
I phase, which is relevant to NaNi02, becomes stable in 
the region where both the Jahn- Teller coupling and the 
Hund'-rulc coupling arc substantial. This phase is un- 
derstood by the strong-coupling picture under the Jahn- 
Teller type electron-phonon couplings. 

3.3 Comparison to experiments 

In our calculations, we successfully reproduce the 
X \/3 COM phase and the t^-OF SF I phase, whose 
ordering patterns arc consistent with the low-T phases in 
AgNi02 and NaNi02, respectively. These two phases ap- 
pear in close parameter regions of 7jt and Jn/U, but for 
different range of the on-site repulsion U. The \/3 x 
COM phase is stabilized in the weakly-correlated region. 
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Fig. 11. (Color online). Energy gap in the r^-OF SF I phase as a 
function of (7 at t' = -1.0, Jh/U = 0.25, 7jt = 1.0, 7br = 1-6, 
and Ai,j. = Ai,j. = 0.05. The line is the simple estimate from the 
strong-coupling analysis. See the text for details. 



where U is smaller than the bare bandwidth, while the 
T2-OF SF I phase is stabilized in the intermediately- to 
strongly-correlated region. 

The difference in U may be attributed to the structure 
of cation bands in these compounds. The magnitude of 
the effective on-site repulsion for Ni 3d electrons is not 
solely determined by its atomic value, but it is consid- 
erably affected by the screening effect brought about by 
electrons in the A cation and oxygen p bands. According 
to the first-principle calculations, Ag bands in AgNi02 
reside in the vicinity of the Fermi level, ^^'^^^ while the Na 
bands in NaNi02 are located about 4eV above the Fermi 
level. Consequently, a larger screening effect is ex- 
pected for AgNi02, which reduces the magnitude of U 
considerably, compared to that for NaNi02. Our results 
are consistent with this trend. In the first-principle cal- 
culations,^^' ^^'^"'^^^ the bandwidth of Ni 'Sd bands was 
roughly estimated to be 2 ~ 3eV, leading to a rough 
estimate of t = 0.25 ~ 0.4eV. In the LSBA+U calcula- 
tion for NaNi02,'^''''^"'^' the value of U was taken to be 
5eV to reproduce the correct size of excitation gap. On 
the other hand, a cluster-model analysis of the photoe- 
mission spectra gave an estimate of J7=7.0eV.'^^-' These 
studies suggest that ?7 ~ 10 — 30t is reasonable, consis- 
tent with our results. It is noteworthy that the gap in 
our calculation at [/ = 20t corresponds to 0.5 ^ 0.8eV, 
which is in the same order of magnitude as the experi- 
mental value ~ 0.24eV in NaNi02.^^ We also note that 
the CO stabilization energy, which is estimated from the 
energy difference between the x \/3 COM phase and 
a para phase, is ~ O.OSi for the present parameters 7br 
and Abr: This result leads a rough estimate that the CO 
stabilization energy is ^ 0.01 - 0.02 eV, which is in the 
same order of magnitude as the CO transition tempera- 
ture observed in AgNi02 (365K26.27))^ 

Finally, we make a brief comment on the absence of 
any explicit ordering in LiNi02. Since LiNi02 is also a 
Mott insulator with a gap of 0.2eV,^"'-' we suppose that 
the compound is in the strongly-correlated region simi- 
lar to NaNi02. In our results, there exists phase compe- 
tition among three insulating phases with different spin 
and orbital patterns, the r^-OF SF, r^^-OF SAF, and Ty- 



OF SF orderings. The competition brings about a frus- 
tration in the spin and orbital sectors in the vicinity of 
the phase boundaries. To argue the consequence of such 
frustration, we need to go beyond the present mean-field 
analysis; however, we can expect severe suppression of 
the orderings and some liquid-like or glassy behavior in 
the spin-orbital coupled system. Hence, one possibility 
is that LiNi02 is located in such competing regime. It is 
noteworthy that the competition is brought about by ex- 
plicitly taking account of the electron-phonon couplings, 
which have not been considered in the previous effective 
model approaches. ■^2' •^^^ In addition to this intrinsic 
phase competition, extrinsic defects on Li sites may play 
an important role in the glassy behavior. Furthermore, a 
long-range strain effect might also play a role through the 
frustrating orbital and lattice sectors. It is interesting 
to take account of these factors explicitly, by extending 
our model and analysis. We leave this problem for a fu- 
ture study. 

4. Summary 

We have investigated the ground state of the multi- 
orbital Hubbard model with electron-phonon couplings 
by the Hartrcc-Fock approximation and the adiabatic ap- 
proximation, in order to elucidate the origin of various 
phases observed in ANi02 in a unified way. We found the 
a/3 X a/3 charge-ordered metallic phase in the weakly- 
correlated region and the orbital-ferro spin-ferro ordered 
insulating phase in the strongly-correlated region. The 
a/3 X a/3 charge-ordered metallic phase is stabilized by 
a compromise between Coulomb repulsions and effective 
attractive interactions from the breathing-type electron- 
phonon coupling as well as the Hund's-rule coupling. The 
electronic state is half metallic; up-spin electrons are lo- 
calized at the charge-rich sites, but down-spin electrons 
are extended and almost free. On the other hand, the 
orbital-ferro spin-ferro ordered insulating phase is stabi- 
lized by the Jahn- Teller coupling under strong electron 
correlation, with a help by the Hund's-rule coupling in 
the competition with other orbital-ordered phases. These 
two phases are promising candidates for the low-T phases 
in AgNi02 and NaNi02, respectively. A possible origin 
of the quite different electron repulsion between AgNi02 
and NaNi02 might be a screening effect from the cation 
and oxygen p bands. The puzzling glassy behavior in 
LiNi02 might be ascribed to the competition among dif- 
ferent spin and orbital ordered states in the strongly- 
correlated region, which occurs under a substantial Jahn- 
Teller type electron-phonon coupling. 
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